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Abstract 

Here we study some general properties of spherical shear-free col- 
lapse. Its general solution when imposing conformal flatness is re- 
obtained [1, 2] and matched to the outgoing Vaidya spacetime. We 
propose a simple model satisfying these conditions and study its phys- 
ical consequences. Special attention deserve, the role played by relax- 
ational processes and the conspicuous link betweeen dissipation and 
density inhomogeneity. 
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1 Introduction 



Gravitational collapse of stars is an important problem of astrophysics and 
building realistic models of collapse remains a formidable task. One inter- 
esting problem is to add heat flow to spherically symmetric models. 

Indeed, dissipation due to the emission of massless particles (photons 
and/or neutrinos) is a characteristic process in the evolution of massive stars. 
In fact, it seems that the only plausible mechanism to carry away the bulk of 
the binding energy of the collapsing star, leading to a neutron star or black 
hole is neutrino emission [3]. 

In the diffusion approximation, it is assumed that the energy flux of 
radiation (as that of thermal conduction) is proportional to the gradient of 
temperature. This assumption is in general very sensible, since the mean free 
path of particles responsible for the propagation of energy in stellar interiors 
is in general very small as compared with the typical length of the object. 
Thus, for a main sequence star as the sun, the mean free path of photons 
at the centre, is of the order of 2 cm. Also, the mean free path of trapped 
neutrinos in compact cores of densities about 10^'^ g. cm. becomes smaller 
than the size of the stellar core [4, 5] . 

Furthermore, the observational data collected from supernovae 1987A 
indicates that the regime of radiation transport prevailing during the emission 
process, is closer to the diffusion approximation than to the streaming out 
limit [6]. 

Many solutions of Einstein's field equations with dissipative fluids carry- 
ing heat flow have been studied (see [7] for references up to 1989 and [8, 10, 9] 
for more recent ones). 

In this vein here we study dissipative spherical collapse with shear-free 
motion. Spherical conformally flat fluids undergoing dissipation in the form 
of radial heat flow where flrst considered in [1] and generalized in [2]. Here 
we reobtain the general conformally flat solution in a slightly different way. 
We match this spacetime to a radiating null held described by the outgoing 
Vaidya spacetime. A simple model is considered satisfying these conditions. 

The paper is organized as follows. In section 2 the held equations are 
presented; in section 3 we reobtain the general solution by considering con- 
formal flatness of spacetime; in section 4 we state the junction conditions to 
the external outgoing Vaidya null radiating held; section 5 presents a simple 
collapsing dissipative model and we flnish with a brief conclusion. 
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2 Field equations 

We assume a sphere of collapsing perfect fluid with heat flow. Its spher- 
ical surface S has center and is filled with radially moving perfect fluid 
conducting heat flow, so having energy momentum tensor 

Tap ^ifJ'+P) WaWp + PQap + QaWp + WaQH, (l) 

where /x and p are the proper density and pressure of the fluid, Wa its unit 
four-velocity, the heat conduction satisfying gato" = and g^is is the 
metric tensor of spacetime. 

We choose comoving coordinates within E and impose shear-free fluid 
motion which allows the metric be written in the form [11] 

ds"" - -A^de + \dr^ + ( dO^ + sin^ ed4>'') 1 , (2) 



where A and B are only functions of r and t. We number the coordinates 



X 



= t, = r, = 6 and x^ = (p and then we have the four-velocity given 
by 

wa = -A5l (3) 

and the heat flows radially, 

= q5f, (4) 

where g is a function of r and t. 

The rate of collapse © = w°'-a of the fluid sphere is given, from (2) and 
(3), by 

e^sA. (5) 

where the dot stands for differentiation with respect to t. 

The spacetime described by (2) has the following non-null components of 
the Weyl tensor Ca/Sjs, 



^4 



^^2323 = y5'sin2^ 



'A' _ B'\ /I B'\ _ fA" _ B"' 



(6) 



and 



C2323 = -r' (|)'sin2^Coioi = 2r^ (^'sin^ ^Co202 

= (I)' Cosoa = -2r'' sin^ OCuu = -2r'Cisis, (7) 



where the primes stand for differentiation with respect to r. 

The non null components of Einstein's field equations Gap = uTa/s, where 
Gap is the Einstein tensor and Tap is given by (1), with metric (2) are 



^2 



2^ _ 
B \B . 



4B' 

2 



(8) 
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B^ 
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-'b-[b 



AB 



npB^ 



(9) 
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~A ^ ^ [^) ^rfi" 

npr^B^ 



B \B 



AB 

^'ab 



Goi — —2 



AB 



A = -kqAB^ 



^From (11) with (5) we obtain 



KqB^ = ^e', 



(10) 

(11) 

(12) 



which shows that the outflow of heat, g > 0, imposes ©' > 0, meaning that, 
if < 0, dissipation diminishes the rate of collapse towards the outer layers 
of matter. If g = then from (12) G' = which means that collapse is 
homogeneous [10]. 

The mass function m(r, t) of Cahill and McVittie [12] is obtained from 
the Riemann tensor component -^23^^ and it is for metric (2) 



m{r, t) = — - — 1X23 



r^B 



p 23 _ 

2 - 2 



B' 



'B^ 



r^B'. 



(13) 



Differentiating m(r, t) with respect to r and t and considering the field equa- 
tions (8-11) we obtain 



m — 



n{rBf{rB)' + gr^B^ 



B 



(14) 
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m — 



K 

2 



pr^B^B + q{rBf{rByA 



(15) 



/^From (14) and (15) we have that the heat flow diminishes the gradient and 
the time derivative of m{r,t). This agrees with the discussion concerning 
(12), since dissipation diminishes the total amount of matter, it is expected 
that the rate of collapse slows down. Furthermore, this agrees too with 
the results obtained for the dynamical instability of nonadiabatical spherical 
collapse [18, 7] where it is proved that relativistically dissipation diminishes 
instability due to the decrease of matter content inside a collapsing sphere. 

The Riemann curvature tensor can be split into the Weyl tensor and 
parts which involve only the Ricci tensor and the curvature scalar. This 
allows to say that the Weyl part is constructed only by the gravitational 
held. Considering the scalar of the Weyl tensor 



with (6), (8) and (13) we obtain after a long calculation 



(16) 



C =48 



m 



XrBY 



KfJ, 



48 



rrir 



{rBf 



where mc is defined as the pure gravitational mass. 



mc^m - '■^iJi{rBf. 



(17) 



(18) 



3 Confer mally flat solution 

Here we impose conformal flatness to the spacetime given by (2), i.e. all its 
Weyl tensor components must be zero valued. Prom (6) and (7) we see that 
if C 22,22, — this condition is fulflUed, hence we have 



'A' 




1 + 2r- 



B'' 



0. 



\A B J \A B J \ B 
We can integrate (19) and after reparametrizing t we obtain 



(19) 



B, 



(20) 
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where Ci is an arbitrary function of t. /,From the isotropy of pressure, (9) 
and (10), equating r^^G22 — Gu to zero and using (20) we find 



which is easily integrated. 



B" B' 1 

2 = 0, 21 

B' B r ' ^ ^ 



^ - City + csity ^^^^ 



where C2 and C3 are arbitrary functions of t. The solution found in [1] is 

a particular case of (20) and (22) with Ci = 0. All conformally fiat perfect 
fluid solutions with g = have been obtained by Stephani [13, 14] 

Conformal flatness imposes C = and from (17) we have that the pure 
gravitational mass mc — and 

m^^f,{rBf, (23) 

which is similar to the result obtained in [15] with q = 0. ^From (14) with 
(23) we have 

fx' = qB'e, (24) 

which shows that for g > and © < then < implying that the density 
diminishes with increasing r. While from (5) with (12) we obtain 

i^f^'^lm', (25) 

which can be integrated, giving 

+ ^^^^ 

where 51 is a function only of t. 

Substituting solution (20) and (22) into (9), (10) and (13) we obtain, 

^^ = ^[%^) ^^^^ 
2(C'2r2 + Cs){C2r^ + C3) - 3{C2r^ + Csf 
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K,p 



{CW + 1)2 



-2- 



Cir2 + 1 



Cir2 + 1 

C2(C2 - 2CiCsy + C^iCCs - 2C2) 



Kq = iiCsCi - C2) 



(28) 
(29) 



Cir2 + 1 

The expansion of the fluid sphere given by (5) with (20), (22) and (29), is 

Kqr Cy + 1 



C2r2 + C3 



4 (C2r2 + C3)^ 



(30) 



We see from (30) that if g = the contraction is homogeneous, however if 
q ^ 0, dissipation produces inhomogeneous collapse, which has already been 
remarked in (12). 

The density /i in (27) confirms the result (26) with g{t) = I2C2C3. 

It is possible to prove, after a long calculation, that the fluid (27-29) does 
not satisfy an equation of state of the form p — cfj,, where c is a constant, 
with q ^ 0. A family of solutions with heat flux satisfying an equation of 
state is given in [16]. 

In the next section we consider the junction conditions of the collapsing 
dissipative fluid to a radiating field. 



4 Junction conditions 



If the collapsing fluid lies within a spherical surface E it must be matched to 
a suitable exterior. Since heat will be leaving the fluid across E the exterior is 
not vacuum, but the outgoing Vaidya spacetime which models the radiation 
and has metric 



1 



2m{v) 



dv'^ - 2dvdp + p^idO"^ + sin^ 



(31) 



where m(u) is the total mass inside E and is a function of the retarded time v. 
In (31) p is a radial coordinate given in a non-comoving frame. The matching 
of these two spacetimes (2) and (31), using the field equations (9-11) and the 
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mass function (13) satisfies [17, 7] 



mivj = < — 



B^B 



[B 



l\2 



B 



r^B' 



iFrom (28), (29) and (33) we have 



3 (Car' + Cs)' Cir\C2r^ + Cs) 



2 C2r2 + C3 



Cir2 + 1 



2(C'3Ci - C2)r 



+2 



(Cir^ + 1) 
C2r2 + Ca 



C2(C2 - 2C,Csy + CsiCCs - 2C2) 



(32) 
(33) 

(34) 



0. (35) 



5 A simple model 

A simple approximate solution for the functions Ci{t), C2{t) and C^lt) sat- 
isfying the junction condition (35) is 



(36) 



where < e <C 1 and a > a constant. When ci = then (36) describes a 
collapsing Priedmann dust sphere, with k — 0, whose radius diminishes from 
arbitrarily large values until, at t = 0, a singularity is formed. The time 
t runs from —00 to and the constant a is proportional to the total mass 
inside the radius r. Substituting (36) into (35) we obtain up to 0(e), 



Cl 



' t 2\ 

Cl f» 0, 



which after integration yields. 



Cl ^ ci(0) exp 



2t 



2ri rs 



Substituting the solution (36,38) into (27-29) we obtain 

12a2 



(1 - e2cir') , 



(37) 



(38) 



(39) 
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1- (1+" '' 



2rv\ 



t Jrl 



(40) 



Kg ^ -e—^, (41) 

which satisfy plausible physical conditions. It should be observed, however, 
that in the general case e 0, the range of t is restricted by physical consid- 
erations. Thus for example if we want the the central pressure not to exceed 
the value of the central energy density, then we should have, 

^>ec,. (42) 

We see from (39) that the energy density diminishes to the outer regions 
due to dissipation; from (40) we have that pressure diminishes too towards 
the outer regions while from (41) we have that the heat flow increases in that 
same direction. 

The mass function (13) inside a radius r with (36) and (38) becomes, 

m(r,t)f» (l-e2cir^), (43) 

showing that dissipation diminishes the mass inside r. Now calculating the 
rate of collapse (5) with (36) and (38) we obtain 

e « ^ (l - ecir^) , (44) 

implying that dissipation slows down collapse. This result agrees with the 
fact that m(r, t) is diminished by dissipation. 
The effective adiabatic index 

r = (45) 

gives a measure of the dynamical instability of the body at given instant of 
time. Calculating (45) for r = and r — r^, with (36) and (38) up to the 
order 0(e) in p we obtain. 



r=0 



2 t 



9 



We see from (46) and (47) that r,_o < ^r=rs which shows that the centre is 
more unstable than the surface region of the collapsing body. This conclusion 
too agrees with our previous analysis. 

5.1 Calculation of the temperature 

Finally it is worth calculating the temperature distribution, T(r, t) , for our 
model, through the Maxwell-Cattaneo heat transport equation [19, 20, 21, 
22, 23], 

rh'^'^w'^qis-j + = -Kh'^^iTp + Tap), (48) 

where r is the relaxation time, K the thermal conductivity and /i"^ = g"'^ + 
w"w^ the projector orthogonal to Considering (2-4) then (48) becomes 

riqBjB + qAB^ = -K(TAy. (49) 

Substituting (20), (22) and (29) into (49) and considering C2 = we 
obtain, up to order e 

kK 1/ 
riCiCsCsJr + C^C^r = -— [T(Cir' + 1)] , (50) 

Now, in the non-dissipativc case ( C\ = eci = 0) it follows at once from (50) 
that T = To(t), implying that in that case the temperature is homogeneous 
within the fluid distribution. Therefore, in the general dissipative case Ci 7^ 
0, we shall have 

T^To{t)+eT,{r,t), (51) 
Then introducing (51) into (50) we obtain up to 0(e) 

rr. rr. / ^\ O 4aVci ( t 2 5\ 9 

where we have assumed for simplicity K =constant and T'c(t) denotes the 
central temperature. The second term on the right hand side of expression 
(52) exhibits the influence of dissipation on the decreasing of temperature 
(remember that f < ) with respect to the non-dissipative cal- 
culated from the non-causal (Landau-Eckart) [24, 25] transport equation, 
whereas the last term describes the contribution of relaxational effects. The 
relevance of such effects have been brought out in recent works (see [16, 26] 
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and references therein). In particular it is worth noticing the increasing of 
the spatial inhomogeneity of temperature produced by the relaxational term, 
an effect which has been established before [27] . 

6 Conclusion 

We have presented the general field equations for a spherical dissipative shear- 
free collapse. Some general properties concerning the effects of dissipation 
on the collapsing body and its mass were discussed. By imposing conformal 
flatness we showed that the system is completely soluble in its radial part and 
producing three arbitrary time functions. Then we matched this solution to 
the outgoing Vaidya radiating spacetime. A simple model with a Priedmann 
limit is constructed satisfying the junction conditions. 

Besides its simplicity, the merit of the model resides in the fact that 
it exhibits in a very clear way the influence of relaxational effects on the 
temperature, and thereby on the evolution of the system. 

It is also worth noticing the appearance of density inhomogeneities di- 
rectly related to dissipation, even though the space-time remains conformally 
flat. This reinforces doubts [28] on the proposal that the Weyl tensor [29] or 
some functions of it [30], could provide a gravitational arrow of time. The 
rationale behind this idea being that tidal forces tend to make the gravitating 
fluid more inhomogeneous as the evolution proceeds, thereby indicating the 
sense of time. 
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